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Abstract Motivated by the diffusion-reaction kinetics on in- 
terstellar dust grains, we study a first-passage problem of 
mortal random walkers in a confined two-dimensional ge- 
ometry. We provide an exact expression for the encounter 
probability of two walkers, which is evaluated in limiting 
cases and checked against extensive kinetic Monte Carlo 
simulations. We analyze the continuum limit which is ap- 
proached very slowly, with corrections that vanish logarith- 
mically with the lattice size. We then examine the influence 
of the shape of the lattice on the first-passage probability, 
where we focus on the aspect ratio dependence: Distorting 
the lattice always reduces the encounter probability of two 
walkers and can exhibit a crossover to the behavior of a 
genuinely one-dimensional random walk. The nature of this 
transition is also explained qualitatively. 

Keywords mortal random walk, reaction-diffusion model, 
first-passage problem, geometry, astrochemistry 



. ^ 1 Introduction 

X. 

^ , Random walks are ubiquitous in statistical physics as well 
- - ' as in the description of many natural phenomena. One key 
example is the study of diffusion-reaction systems. 

In our particular case, one motivation for the following 
examination comes from a long-standing astrophysical puz- 
zle, namely why there can be a large abundance of molecu- 
lar hydrogen in interstellar clouds of gas and dust compared 
to the concentration of atomic hydrogen. For realistic condi- 
tions, the generally accepted mechanism is as follows [1^]: 
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Reactants impinge onto the surface of dust grains in a homo- 
geneous fashion. They diffuse on the surface. Additionally, 
they desorb at a certain rate. If two such atoms meet during 
their traversal of the surface, they react to form a hydrogen 
molecule which immediately desorbs. We will sometimes re- 
fer to the recombination efficiency of such a system, defined 
as the fraction of impinging atoms that eventually react to 
yield a molecule. 

To study the catalytic role of the dust grains in the as- 
trophysically relevant regime of a small average number of 
atoms on the grain, one sets up a zero-dimensional master 
equation for the distribution of the number of reactants on a 
single grain, including terms for deposition, desorption and 
reaction, and solves for the stationary state of the system 
[5,6]. Rate equations for the mean adatom number do not 
suffice as they cannot account for crucial fluctuations and 
the non-PoiSSONian statistics of the reactant number that 
is induced by their confinement to a finite closed surface 
on which the reaction takes place [6,7]. Within such a mas- 
ter equation treatment, apart from obvious single-atom rates, 
there only occurs a single two-particle parameter, viz. the 
sweeping rate describing how often a pair of atoms on the 
grain meets. One convenient way to express this parameter is 
by means of the encounter probability originally introduced 
in [8] in a continuum setting. This is the key quantity of a 
first-passage problem of two mortal random walkers: Mor- 
tality (i.e., terminating a random walk with a certain proba- 
bility per step) corresponds to desorption of the atoms on the 
surface, the random walk to their diffusion, and the question 
(of first-passage nature) is whether the two walkers meet be- 
fore either one dies, quantified by the encounter probability. 

In general, a realistic full diffusion-reaction model has to 
be examined by microscopic Monte Carlo simulations (see 
e.g. [9]). This is mostly owed to the fact that applications, 
and especially the astrophysical problem we described, call 
for the inclusion of disorder in the local rates of hopping 
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and desorption, a feature hard to tackle analytically in more 
than one dimension. Despite the vast literature on diffusion 
properties in disordered media this even holds true (to the 
best of our knowledge) when we focus on the basic level of 
the encounter probabihty. 



Still, it is useful to start with the analysis of simple ho- 
mogeneous models that are accessible rather easily, both to 
establish a reference point for the validity of simulations and 
to make progress in the analytical exploration. On the level 
of the encounter probabihty this means we can employ re- 
sults from the theory of homogeneous random walks. Fur- 
thermore, the discrete random walk model allows us to ex- 
tend the generality step by step to different lattice types and 
shapes, hopefully providing deeper insight into the impor- 
tance of these factors. Finally, it is expected that the problem 
will have some fundamental appeal for the theory of random 
walks in relation to reaction-diffusion systems in general. 



Examples from seemingly remote fields to which this 
analysis might bear relevance include second-layer nucle- 
ation in epitaxial crystal growth [8], chemical kinetics in- 
side aerosol droplets [10], biophysical problems like exciton 
trapping on photosynthetic units [1 IJ and a range of search, 
transport and binding processes of or along DNA strands 
[12, 13]. Admittedly, in the latter case the homogeneous sit- 
uation is only a first step to meaningful results, and the cru- 
cial two-dimensionality of our problem does not necessarily 
apply either. 



Our goal in this paper is thus to further elucidate the 
meaning of the encounter probability, examining the pitfalls 
of the continuum limit to aid comparison to Monte Carlo 
simulations, and to analyze the influence of the lattice type 
and, most importantly, different geometries. We strive for a 
completely analytic theory to fully explain all findings of 
simulations. To this end, we will proceed as follows: Sec- 
tion 2 will first introduce the model and basic definitions 
along with our notation. We will then show how the funda- 
mental quantity that we call 'encounter probability' can be 
obtained from a simple random walk calculation, and will 
analyze its asymptotic behavior. In Section 3 we investigate 
the subtle continuum limit which implies a logarithmically 
slow convergence to results obtained in a continuum model, 
which is crucial for comparison to simulations. We supple- 
ment this by simulation results in Section 4. Section 5 we ex- 
tend the random walk calculations to a rectangular lattice to 
examine the role of the lattice shape, analyzing a quasi-one- 
dimensional limit and deriving and explaining the effect of 
a distorted aspect ratio. FinaUy, we present our conclusions. 



2 Random walk treatment 

2.1 Model and definitions 

We consider a homogeneous two-dimensional lattice with 
periodic boundary conditions, meaning that all S = Li Li 
sites are equivalent. Two walkers are randomly deposited on 
two (not necessarily different) sites of the lattice, then they 
move between nearest-neighbor sites with an undirected hop- 
ping rate a, and may desorb at a constant rate W <^a. When- 
ever the two meet by one hopping onto the site occupied 
by the other (or by coinciding initial positions), they react 
and this random walk reaUzation has ended successfully; the 
probability (averaged over random walk realizations and ini- 
tial positions) at which this happens before either of the two 
walkers desorbs is called the encounter probability p (sub- 
scripts will denote certain models and /or simulations from 
which it is obtained, as well as further specializations). 

Here we will treat a discrete- time random walk, although 
a continuous-time version (CTRW) with an appropriate wait- 
ing time distribution might seem closer to the physical sys- 
tem at hand. We do this in order to keep the analytical treat- 
ment as simple as possible, and we will later argue and nu- 
merically prove (Section 4) that the discrete version accu- 
rately describes our situation. Due to homogeneity of the 
lattice, the problem in fact reduces to that of a single walker 
meeting a target site, cf. Appendix A. 

2.2 Exact results 

Quite generally, we are concerned with finite lattices transla- 
tionally invariant in all directions, which extend to Lj lattice 
sites in the Jth of d dimensions before periodically contin- 
uing. The exact expression of the random-walk encounter 
probability on such a lattice is re-derived in Appendix A and 
reads [14] 



p^ = {\-^)SP*{Q;^) 



(1) 



mei3 



^A(2a:L-im)' 



Here, 



1 



<1 



(2) 



is the survival probability per step, S = nf=i Lj is the to- 
tal number of sites, and P* (0; ) is the number of times a 
mortal walker on a periodic homogeneous lattice returns to 
the origin. Q denotes the lattice, m is a 'lattice vector' of 
d integer components < mj < Lj — 1 with j = l,...,d, 
and L = diag(Li .... ,Lj). Finally, A is the structure func- 
tion (basically a discrete Fourier transform of the normal- 
ized transition probabiUty) of the walk. We specify to J = 2 
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at this point, where it reads A(k) = (cosfci +cosfc2)/2 for 
an isotropic walk on a square lattice that we will further on 
label as being of 'type (a)', and A(k) = [cos^i +cosA:2 + 
cos(^i +^2)]/3 for the isotropic walk on a triangular lattice 
(coordination number 6), now designated as 'type (b)'. 

We focus on the two-dimensional case with both lattice 
lengths much larger than unity and with 'long survival' de- 
fined by 1 — ^ 1. The expression (1) then affords several 
regimes, characterized by the comparison of dimensionless 
'lengths' . Introducing the typical single-atom random walk 
length 



1-^ 



> 1 



(3) 



and with lattice dimensions Li 2 ^ 1, one can associate 1 <C 
£ <C L12 with 'large' lattices, and 1 <C ii.2 "C i with 'small' 
lattices. The intermediate regime in which one lattice length 
is smaller, yet the other larger than the random walk length 
is to be discussed in a later Section. 

For a still fairly large class of walks that includes the 
two cases (a) and (b), one summation in equation (1) can 
be carried out explicitly [11]. Appendix B gives the result, 
which we generalized to include the case Li ^L^- Whenever 
we need to numerically evaluate prw we use the resulting 
single-sum expression, further simplified for the two lattice 
types, and implemented in a small GNU Octave script. 

The encounter probability as given above allows for the 
two random walkers to start on one site, and this counts as 
an encounter on the zeroth step. For applications and com- 
parison to other models or simulations, we will often need to 
deal with the encounter probability calculated such that 
it only accounts for meeting of the walkers by hopping, and 
that does not allow the initial condition of walkers starting 
at the same site. The two quantities are related by 



/'rw = 7 + ( 1-7 JPrw, 



(4) 



as becomes obvious by splitting up according to the 
starting site of the second walker, i.e., either on the same 
site as the first one (with probability 1 or on any other 
site (with probability 1 — 1/5). This leads to the expression 
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in terms of p™, and p with other subscripts will henceforth 
denote probabilities p that are obtained from other models 
using the corresponding analogous convention, i.e., exclud- 
ing an encounter due to coinciding initial positions. In the 
context of hydrogen recombination on dust grain surfaces, 
the corresponding mechanism is the 'Langmuir-Hinshelwood 
rejection' of atoms that impinge on top of another; see e. g. 
[15] and references therein for early original work. 



2.3 Large lattice approximation 

Large lattices, formally given by 5 » a/W » 1 or, more 
precisely, by 1 <C ^ <C ii,2, can equally well be defined as 
the regime in which boundary conditions, and (apart from 
the initial placement) the overall number of sites of the lat- 
tice no longer matter at all - even if the boundaries truly 
affected the random walker (say, by reflection), it either ap- 
proaches the target or experiences the finiteness and (possi- 
bly) boundedness of the lattice, but never during one single 
walk. It is thus safe to simply send Lj 2 — > 0° in f*(0;i^) 
in (1), by which procedure the sum becomes an area integral 
(/"(O;^) in the standard notation used in the Appendices). 
Some more details on the ensuing approximations are given 
in Appendix C, with the result that 
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square lattice, 
triangular lattice. 



(6) 



with relative errors of ^(1 — 1^). 



2.4 Small lattice approximation 



Here we have 1 < Li 2 <C ^ or a/W > 5 > 1. Again, we 
employ results for P* (0; ^ ) from the random walk literature. 

For the moment, we restrict ourselves to the case L\ = 
L2 = Vs. The expansion for 1 — § ^ 5^^ ^ 1 then reads 



1 



+ ci]n{cS) 



(7) 



where c and ci are real constants. This is shown in [11], also 
presenting the first calculations that deliver upon the impor- 
tant pre-factor inside the logarithm (for square and triangular 
lattices), which were extended and subject to minor correc- 
tions in [14] (cf. the earlier and easily accessible derivation 
in [16], which unfortunately does not give this pre-factor). 
For the encounter probability one thus obtains 



SW 

Prw ~ 1 ciln(c5). 



(8) 



It is important to note that the original expansion is valid 
in the regime (1 — i^)57r~'ln(c5') <C 1, or equivalently, 1 — 
Prw ^ 1 . This would thus be a more precise definition of a 
'small lattice'. 

The constant ci has the value I/tt for the square lat- 
tice case (a), and ci = \/3/ (27r) for the triangular lattice 
(b), respectively. The crucial factor c inside the logarithm 
appears in the different guise C2 in [11], related to ours by 
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cijcx = Inc. For the square lattice the ratio yields 



n 1 
lnc= -+2(7-ln7r + -ln2) 
3 2 ' 



+ 4 

w 0.612807020, 
while for the triangular lattice this becomes 



Inc : 



^ + 2(r-hi;r + ihi3) 



e 2 ^3 



! 0.853262084, 



(9) 



(10) 



where we used standard rounding and checked the consis- 
tency of the numerical evaluation against the original figures 
and, when available, improved ones from [14]. Here and in 
what follows, 7 w 0.5772156649 denotes Euler's constant. 



We will also refer to results for p obtained from random 
walk theory, however in a heuristic fashion [7], using known 
expressions for N^is, the average number of distinct sites vis- 
ited by a random walker after it has taken a given number of 
steps. The asymptotic result is that after I steps (on a 
two-dimensional regular lattice), 

nn 



(11) 



' C\n{Bn) 

with real positive constants C and B depending on the lattice 
type, viz. C = 1 and 5 = 8 for the square, and C = a/3/2 and 
B = 12 for a triangular lattice, respectively (see e.g. [17]). 
The encounter probability derived on this route will be la- 
beled Prw.hcur- Note that the above A^dis is the leading term 
of an asymptotic expansion, and thence a poor estimate for 
moderate values of n. 

In the comparison of analytical results, we allow for walk- 
ers to meet by coinciding initial positions throughout, and 
we incorporate corresponding terms in the diffusion models 
as well; in short, we always use p and not p. 



3 Continuum limit 



3.1 Large lattices 



A natural question to ask is whether the model described 
above yields a reasonable continuum Umit for Li 2 — » °°; on 
the one hand, since such continuum limits are of genuine 
theoretical interest in themselves, on the other hand, since 
we have earUer solved a continuum model [7] to compare 
against. The answer will also provide important information 
for the comparison between analytical theory and simula- 
tions. 

The proper scaling of the discrete random walk parame- 
ters then proceeds as follows. We want to stay in the same 
regime of the system, which means we need to keep both ra- 
tios Lj/£, and consequently SW/a, fixed. Further, we do not 
want to distort the given aspect ratio of the lattice, so that 
Li /L2 is kept constant as well (a moot point as long as we 
use quadratic lattices, i.e., Li = L2 anyhow). With these con- 
straints we let the number of adsorption sites become very 
large, S ^ °°. Such a joint hmit preserves all quantities in 
our expressions that only depend on the 'regime parameter' 
~ 5(1 — ). The limiting behavior can thence be determined 
from the corresponding approximate results for the two re- 
gimes separately. 

In the continuum (or "diffusion") model proposed in [7], 
a stationary diffusion equation for the probability density of 
a moving atom is solved on a spherical surface, and given 
an absorbing 'target area' . The encounter probabiUty is then 
calculated as the (suitably normalized) diffusion current en- 
tering this area. The results will be denoted by pdie-^ 

' Note that the cited article differs in that it uses this notation for 
only a certain contribution to the full encounter probability. 



The random walk result pnv has been given above. Its diffu- 
sion model analogue reads 

4ga 1 



PdiS' 



SW ln(4fl/W)-27' 



(12) 



where g denotes a lattice-dependent factor; it is g = 7r/4 for 
the square lattice, and g = K/ (2\/3) for the triangular lattice, 
respectively. The heuristic result obtained from NcHs is 



Prw,heur ' 



na 



1 



CSW ln{Ba/W) 



(13) 



with B, C as defined above. Since comparison between the 
factors g and C immediately shows that for our cases 4g = 
n/C, we see that not only the functional dependence, but 
also the numerical pre-factors of all three expressions prw. 
PdiS, and /Jrw,heur Coincide. 

As far as the factor inside the logarithm is concerned, it 
is not too surprising that pi^, and /Jrw.heur agree as well. How- 
ever, the diffusion model expression pdis clearly differs here, 
though this becomes irrelevant in the true limit 5 — > t». Note 
that there is no longer any imprint of the lattice type, as op- 
posed to the discrete result. We should emphasize that in [7] 
care was taken that in the given form of the asymptotic ex- 
pressions, all omitted terms are of higher order, particularly 
in the denominator expression, where further terms are of 
'polynomially' smaller order than unity. 

With this, the relative error of p^iff with respect to the 
random walk expressions can be calculated to yield 



A = 



lnB-2(ln2-7) 
hi(a/W)+2(ln2-7) 



(14) 
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Since 2 (In 2 — 7) ~ 0.23 and taking into account the values 
of B this means that in the large lattice regime, the diffu- 
sion model encounter probability (and in its trail the sweep- 
ing rate and the recombination efficiency) systematically ex- 
ceeds that given by the random walk result. The discrepancy 
vanishes in the continuum limit 5 — > t», however it only does 
so like \/\n{a/W) or as \/\nS (since SW /a = const.). 



3.2 Small lattices 

The earlier results to compare the prw asymptotics against 
are as follows: For small lattices, we obtained the diffusion 
model expression 



SW 

4ga 



(15) 



with g as defined in the previous Section. The heuristic ran- 
dom walk derivation yielded 



CS^V 

Prw.heur « 1 ln[BC/K ■ S] . 

na 



(16) 



Including the factor c\ occurring in p^sN of (8) in the relation 
established between g and C, we observe that it satisfies 



1 



C 



\g n 



(17) 



Thus once again, the pre-factors as well as the functional 
dependence of all three expressions agree. 

Turning to the pre-factor inside the logarithm, we have to 
resort to the numerical values of Inc as given in Section 2.4. 
Starting with the /^rw.heur expression, numerical evaluation of 
the corresponding XnBCjn yields numbers of 0.934711656 
and 1.196335728 for type (a) and (b) lattices, respectively. 
While not in good agreement with the values for Inc, this 
is still reasonable: Apart from the nature of the A'dis expan- 
sion, the ensuing heuristic derivation asserts that the walker 
does not die until it has traversed all S sites, and includes 
inversion of the transcendental equation Ndis{t) = S which 
is approximated by one iteration. It is the very pre-factor 
just mentioned that suffers from this approximation. In con- 
trast, the analogous factor in the diffusion model result reads 
— In^— 1, and numerical evaluation provides —0.758435525 . 
and —0.902276561... for types (a) and (b), respectively, 
both completely off. 

For S °° with SW /a = const, we will eventually still 
leave the small lattice regime as p drifts away from unity, cf. 
our earlier remark regarding the used expansion. To measure 
the error of the diffusion model, it is hence more appropri- 
ate to compare the complements I — p of the two models, 
leading to a relative deviation 



A' = - 



In(cge) 
ln(c5) 



<0 



(18) 



of the diffusion model outcome relative to that of the ran- 
dom walk. Obviously, again the diffusion model deUvers 
an encounter probability systematically larger than obtained 
from the random walk model, and again, this discrepancy 
vanishes only logarithmically with the system size S ^00. 



3.3 Results 

Besides the interest in the exact encounter probability of 
the random walk p^w and its limiting behavior for the two 
regimes, we now have estabUshed that its asymptotics dif- 
fer in logarithmic terms from that of the continuum model 
Pdis- While the difference between discrete model results 
and those of the continuum model vanishes in the true con- 
tinuum limit S ^ °°,W /a ^ °° with SW/a = const., it does 
so only logarithmically in the system size S: This slow con- 
vergence is an inevitable direct consequence of the marginal- 
ity of spatial dimension two for the random walk and diffu- 
sion, unfortunate for the comparison of discrete-space simu- 
lations to analytic results. 

We finally plot (Figures 1 and 2) the relative discrepancy 
of the exact results of the diffusion model pdiff with respect 
to that of the random walk />rw per cent, i.e., 100 • (pdis — 
Ptv/)/Piv/ for the large lattice regime, substituted by the rel- 
ative difference of the complements (1 — p) for the small 
lattice plots. Inside a single plot, we sample lattice sizes de- 
termined by the lengths Lj = L2 = L being the closest inte- 
gers to give 5 = 4 • 102, 4 • 10^ . . , viz. L = 20, 63, 200, 632, 
2000, 6325. Simultaneously, \og{a/W) =2,3... increases 
with lattice size, such that the 'regime' parameter SW/{4a) 
stays nearly constant and we properly approach the contin- 
uum Umit. To compare, we also plot the error estimates A 
and A' as appropriate for the given regime. 

In both regimes the error predictions nicely agree with 
the discrepancy of the exact results, with increasing preci- 
sion the further we go into a certain regime. This also veri- 
fies that the asymptotic results used to derive the estimates 
are correct, and that the discrepancy is not due to any rele- 
vant terms erroneously omitted. The error generally shows 
the predicted slow logarithmic decrease in the system size. 
Moreover, one can clearly see that it is considerable for all 
lattices, even for the largest ones, and the error even increases 
the further we enter one regime. 



4 Simulations 

We already mentioned that in many applications, the full mi- 
croscopic kinetics is simulated within the Monte Carlo ap- 
proach. Thence here we quantitatively compare our analytic 
predictions for the encounter probability to the correspond- 
ing simulation results. It is fairly easy to simulate both, the 
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Fig. 1 Relative difference (in percent) of 1 — p^iff w.r.t. 1 — p^x for 
SW/(4a) = 10"^ (top) and SW/(4a} = IQ-^ (bottom, still a imaH lat- 
tice), type (a) (dashed / squares), type (b) (dotted /triangles); thinner 
lines the prediction — A', 
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Fig. 2 Relative difference (in percent) of pas w.r.t. prw for SW/(4a) = 
10° (top, already a large lattice) and 5W/(4a) = 10^ (bottom), type (a) 
(dashed /squares), type (b) (dotted /triangles); thinner lines the predic- 
tion 4. 



stochastic two-particle system for which we defined the en- 
counter probability p, but also the 'full' diffusion-reaction 
system including continuous stochastic deposition of new 
random walkers, yielding results for the recombination ef- 
ficiency of this process. 

We use a standard kinetic Monte Carlo algorithm to keep 
track of individual atoms deposited onto, hopping around on 
and desorbing from the lattice. As in the random walk de- 
scription, the lattice is homogeneous, with periodic bound- 
aries. All waiting times for events are stochastically distributed 
according to an exponential distribution with the average 
waiting time given by the inverse rate of the corresponding 
process. In the simplified case of the two-particle simulation 
we consider two independently moving and desorbing atoms 
in continuous time - though with an arbitrarily rescaled time 
step permitted by the homogeneous envirormient. This is 
the simplest implementation of the continuous-time random 
walk (CTRW) introduced by [16]. 

The algorithm respects the Langmuir-Hinshelwood re- 
jection described earUer, i.e., atoms impinging on top of 
occupied sites are repelled. Hence for faithful comparison 
between simulation and random walk theory, in this whole 
Section we will refer to prw instead of p^vi- 

Let us first relate what follows to earlier work. The sim- 
ulation of a single moving (and possibly desorbing) atom 
that is to meet a fixed immortal target site in a homogeneous 
environment (hence the situation used in deriving prw) obvi- 
ously obliterates any notion of continuous time and can be 
fully described in terms of steps of the walk. It is therefore a 
truism that any correctly working microscopic Monte Carlo 
simulation of such a system has to converge to p^^ if we sam- 
ple enough individual trials. For the case without desorption 
(where the interest is not in the probability of an encounter, 
but rather in the statistics of the number of steps this takes), 
the predictions of [1 1] have been confirmed by Monte Carlo 
simulations in [18], where the effect of alternative boundary 
conditions of the finite lattice is also examined. 

Now if both reactants can move and desorb, it makes 
sense to simulate them in continuous time to let all events 
happen in an ordered fashion, and this is what we have done 
in all simulations we refer to herein. The equivalence of 
both scenarios (two walkers vs. one walker, and regarding 
the encounter probability) was clear in the stationary dif- 
fusion model, due to its continuum nature, the absence of 
any time, and its dealing with a "concentration" that is ob- 
tained by averaging over random walk realizations. For the 
single instance of the CTRW we note that (in our homoge- 
neous setting) we can still treat one of the two walkers as 
fixed: If it is to make a move, we may simply re-label lattice 
sites accordingly and let the other walker perform an appro- 
priate step instead. If it dies, the trial is ended as well as 
if the other walker had died. Hence the system can equally 
be described as a single walker that is to meet a target site. 
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hopping and dying at twice the rates of each of the two orig- 
inal random walks, and still with exponentially distributed 
waiting times. This change of rates does not alter the sur- 
vival probability per step and thence does not change the 
encounter probability. Moreover, while the number of steps 
performed in a given time is a stochastic quantity for the 
CTRW (as opposed to the discrete-time random walk, their 
relation for a huge class of waiting time distributions exam- 
ined in e. g. [19]), the time passed never matters for the en- 
counter probabiUty, so that a description in terms of steps 
only is fully justified. 

Indeed, we can report excellent agreement of the encounter 
probability prw obtained from the random walk model with 
that found in Monte Carlo simulations, pmc, see Appendix D 
for details. This holds for arbitrary parameters S and W/a, 
as well as for both types of lattice, and gives credence to the 
correctness of our simulation method. 

To obtain p^c we sampled A' =10^ individual trials (ex- 
cept for the largest S and £, where we chose = 10^ due 
to time constraints) of the Monte Carlo simulation of two 
atoms moving in continuous time. Such a repeated sampling 
of independent Bernoulli trials with (unknown) exact suc- 
cess probability p (corresponding to recombination in our 
context) yields an absolute standard deviation of the out- 
come of a = \/p{l —p)/N. Thus a = a/p, evaluated using 
the (known) values for prw> is the expected relative standard 
deviation for simulations. Obviously a becomes fairly large 
for small p, i.e., approaching the large lattice regime. 

The binomial distribution for the number of recombina- 
tions with fixed success probability p of the individual trial 
and for the number of samples N ^ °° tends to a Gaussian 
distribution, so the distribution for the numerically found 
pmc should tend to a normal distribution with the above stan- 
dard deviation. Roughly a third of all results are then ex- 
pected to he outside a corridor of half width a due to fluctua- 
tions; counted over all simulation data points for the quadratic 
lattices we find precisely 26/78 values outside this range. 
The largest deviations from p^^ are of the order of less than 
3a, with the largest absolute deviations smaller than 7%. 
Furthermore there is no systematic over- or underestimate of 
the theory. Similar results were obtained for rectangular lat- 
tices (see Appendix D). We conclude that the simple random 
walk result describes the continuous-time simulations to ex- 
cellent accuracy. On a side note, this shows that the diffusion 
model results would completely fail to describe simulation 
results, most pronounced in the large lattice regime. 

The excellent agreement between simulations and ran- 
dom walk theory allows us to revisit the effect of the lat- 
tice type on the recombination efficiency, which was dis- 
cussed in [7] with reference to the Monte Carlo simulations 
reported in [9]. We have seen above in Section 3 that, to 
leading order, the lattice type enters the encounter probabil- 
ity through the multiplicative constant (17), which differs by 
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Fig. 3 Random walk encounter probability p„, as a function of 
the (logarithmic) grain size for W/a = 10^', and for a type (a) 
(dashed /squares) or a type (b) (dotted /triangles) lattice, respectively. 



about 15% between the square and triangular lattices. The 
effect obtained using the full, exact expression for is of 
the same order or smaller, as illustrated in Figure 3. Since 
neither the sweeping rate A nor the recombination efficiency 
depend more strongly on the lattice type than p itself [7 J, we 
may conclude that significant lattice type effects cannot be 
expected in simulations that faithfully represent the reaction- 
diffusion model treated in this paper. The discrepancy with 
the results of [9], where efficiencies on square and triangular 
lattices were found to differ by a factor of 2, thus remains. 



5 Shape of the lattice 

As was mentioned in [7,8], the diffusion model results for 
the encounter probability can be compared to those obtained 
in a completely analogous way for a flat disc with a reflect- 
ing outer boundary and a fixed absorbing target in its cen- 
ter. A detailed analysis (again also accounting for encoun- 
ters "by deposition") shows that upon natural identification 
of the parameters, in both regimes the functional form and 
pre-factors coincide, but while in the "large-disc" regime 
the asymptotics reproduces the sphere result (12) exactly in- 
cluding the logarithmic factor (or coequal ^(1) terms), it 
slightly differs in the "small-disc" regime, where 

SW 

Pdm^isc ~ 1 - HS/g) - 3/2] , (19) 
4ga 

with the lattice factor g defined before. Numerical evalua- 
tion shows that there is only a small (of the order of a few 
percent) difference of both models in between these hmiting 
cases. 

This re-assures our conviction that curvature effects should 
not matter with S ^ I: The basic reason is that the radius 
of curvature (in units of lattice spacing) is of the order of 
the system size. The absence of any effect is obvious for 
the large-lattice regime, when a random walk typically does 
not travel a distance long enough to feel the radius of cur- 
vature, £ <C Vs. But as soon as £ ~ v^, nearly the whole 
lattice is swept by the walk anyway, the smaU fraction not 
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explored not depending on curvature, but on the failure of 
locally dense exploration. 

Note that the difference of the logarithmic factors of 
sphere and disc in the small-lattice regime appears in the 
same order of magnitude as the discrepancy between lattice 
and continuum models in Section 3, although the former 
models (in contrast to the latter ones) substantially differ 
in topology and boundary conditions. While at some point 
we had considered curvature and connectivity effects respon- 
sible for the discrepancy between the diffusion model re- 
sult and that obtained for the random walk, this comparison 
strongly (and quantitatively) suggested otherwise. 

It is clear, however, that boundary conditions and the 
general shape of the lattice should have some effect on the 
random walk properties (at least in certain parameter regimes), 
and possibly on the encounter probabihty p. Related ques- 
tions concerning the importance of confined geometries, with- 
out desorption and correspondingly focusing on mean first- 
passage times instead of first-passage probabilities, have been 
examined in detail in a recent remarkable series of papers 
[20-22]. Moreover, as we have seen, the peculiarities of spa- 
tial dimension two are very important for the features of p 
- what happens then, is one naturally ensuing question, if 
we distort the lattice shape in such a way that it becomes 
effectively one-dimensional? 

In view of the line of thought of this article we focus here 
on the analysis of a torus lattice (periodic and rectangular) 
with distorted aspect ratio, that is, with different lengths Li 
and L2 in the two directions. We will refer to the case Li ^ 
L2 as rectangular, while we call Li = L2 a quadratic lattice, 
the latter not to be mixed up with the term square that we 
use exclusively to label a certain lattice type (i.e., its internal 
structure) as opposed to its shape or geometry. 

Let us therefore refine our notion of small and large lat- 
tice regimes. We still assume all three lengths Li, L2, » 1 
throughout. Without any loss of generality, let L2 < Z-i, but 
we might further specialize to L2 *C Li . Ordering lengths, we 
are then left with three (instead of our former two) asymp- 
totic regimes defined by the ordering of length scales: 

• K£<L2 <i'i. 

• 1 «:z,2<^«:Li. 

• 1 «: L2 < i-i «: ^• 

The first case corresponds to the earlier large-lattice regime; 
the third case to small lattices. The intermediate regime in 
which the walk easily sweeps one dimension but is typically 
short compared to the other is a new property. We define the 
aspect ratio as /i := L-ijLx G (0, 1]. 

We will dismiss the large-lattice case for most of the se- 
quel. In this regime, boundaries are effectively not felt by 
the random walker, and the appropriate double integral limit 
of (1) yields a result that only depends on the total number 
of sites S (and no longer on Li^2), as given in Section 2.3. 
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Fig. 4 Encounter probability (on a square lattice) as a function of log ^ , 
for fixed sizes 5 = 4x10^ (thick) and 5 = 4x10^ (thin), respectively, 
mASW l{^a) = 10"^ in both cases. 
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Fig. 5 As Figure 4 with SW/(4a) = 1. 

5.1 Overall behavior 

We keep 5 and ^ constant and show the encounter proba- 
bility prw on a type (a) square lattice as a function of (the 
logarithm of) the aspect ratio [i for both, one absolutely 
small (5 = 4 X 10^) and one absolutely large lattice (S = 
4 X 10^). In Figure 4, we start in the small-lattice regime for 
the quadratic case (rightmost in the plot), whereas Figure 5 
starts in the large-lattice regime. There is no qualitative dif- 
ference between the two lattice types with respect to any 
plot. 

In all cases the encounter probabihty shows a strong 
monotonic decline upon distortion (i.e., moving left in the 
Figures), though with a less pronounced shape and ending 
at a fair fraction of its peak value for the absolutely small 
lattice. We will first analyze several aspects of this transition, 
finally putting our findings together in a qualitative explana- 
tion. 



5.2 Extreme distortion 

We start with a simple observation. First, let us use the two- 
dimensional result for p^v on a type (a) square lattice, and 
let one lattice dimension shrink to one lattice unit only (S = 
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LiL2 =L \). What we get from (1) is 



m=0 



(,^/2)(cos(27rOT/5) + l) 



(20) 



where we have not yet used any summation results. Second, 
we immediately start with a genuine one-dimensional lattice, 
and then the sum for an unbiased walk (different structure 
function!) reads 

I I . J--,.. , (21) 

A ^ 



^'rw'.ld 



t, co&{2nm/S) ' 



which is obviously different from the former hmit. We see 
that prw, Id > Piw for equal number of adsorption sites S and 
equal ^ . 

The explanation of this result of an "extremely distorted" 
situation provides some insight for the general problem. In 
the two-dimensional version, the walker still performs steps 
in the direction in which there is only one lattice unit length 
(of the periodic lattice), and one can imagine it as walking 
around on a very 'thin' torus with one dimension completely 
"wrapped up". Contrary to this, the truly one-dimensional 
walk only performs steps in the direction in which the lat- 
tice is extended. The net effect is that the two-dimensional 
walker wastes on average half the number of steps it takes by 
walking in the 'wrong' dimension and not coming any closer 
to the target site. In fact, this suffices to re-gain one result 
from the other heuristically: The 2d-walk; corresponds to a 
Id-walk with the same desorption rate W, but with the undi- 
rected hopping rate effectively halved by the useless wait- 
ing steps, a — > a/2. This change translates to (1 — ^)/^ — > 
2(1 — ^)/^, and minimal manipulation of the two expres- 
sions given above shows that this casts the Id- walk result 
into that for the two-dimensional walk (see [23] for a related 
discussion). 

5.3 Small lattices and Li Li 

For the following analysis, we generalized the asymptotic 
resuh for P* (0; ) of [1 1, 14] to the case Li ^ Lj, at least in 
the vicinity of the square lattice. To this end, one mainly has 
to diUgently separate the different L^s in the derivation and 
to check that the individual steps continue to hold. Based on 
the generalized exact expressions given in Appendix B, we 
obtain 

1 InLi 



P*(0;<^): 



LiL2(l-<?) '-;r(l-2^o) 



W + ^ [2y+ 21n(2/;r) - ln( 1 + t] )] + ^"Vr 
2(1 -2^o) 

(3i)-l)jiL, , ^(1) 



(0), 



+ 



+ 



(22) 



36r 



2(1 -2^o)LiL2 

+ff{L-^) + ff{\-^yi^, 



where, first of all, terms are arranged in orders of ^J\ — ^, 
as this is the smallest quantity, and, inside, according to or- 
ders of Lt^2- The parameter is the probability a step is 
directed into a particular lattice unit direction, qi and q2 de- 
note transition probabihties into the diagonal directions, cf. 
Appendix B. Furthermore, we adopted the definitions 



J] 



2[(go + 2gi)(go + 2g2)]^/' 
1-2^0 

go(l -2go) 
{qo + 2qi){qo+2q2) 



(23) 



1. 



which take the values r = 1 = tj for the square lattice case 
(a), and r = rj = 1/3 for the triangular lattice case (b). 

The only parts in the above expression that have to be expUc- 
itly evaluated for the two lattice types are the sum contribu- 
tions sf'^\ Unfortunately, their generalization to Li ^ L2 is 
reasonably feasible only for the square lattice case (a), on 
which we now focus and where 



^3 - nV +2^ +3 

= (^^-2KI^/L, +3e-4;rL2/L, 

+12^'^"^^^^ +...y 



(24) 



Due to the nature of the involved approximations, the 
expansion is asymmetric in the lattice lengths. In the form 
given above, L2 is the Umit of the 'inner' sum that is ex- 
plicitly and exactly evaluated (see Appendix B), while the 
other length L\ appears as the upper limit in several Euler- 
MacLaurin formulae in the derivation. Therefore we expected 
the result would fit best for the case where L2 < Li . This is 
the ordering we will consistently assume throughout the fol- 
lowing. 

The detailed assumptions on the length scales that enter 
the result are as follows: First, we need LiL2(l — ^) ^ 1 
for the general form of the P* expansion to be meaningful, 
cf. [11], Eq. B35. Further, the result rehes on the fact that 
the dimension for which the sum is explicitly evaluated is 
small, L2<^i= — « 1 / \/\ — t„ necessary for the 

used expansion of the first summand to converge {ibid. Eqs. 
Bll, B12), while the other is absolutely large Li ^ 1 (for 
applicabihty of the Euler-MacLaurin formula). Lasdy, 
L2/L1 must not be too small (owed to convergence in Eqs. 
B28 ff. ibid., also see (24)). This set of conditions appUes to 
the small-lattice regime whenever L2/L1 does not become 
too small; but it is also satisfied for the intermediate regime 
provided that additionally, LiL2(l — 1)^1- 
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It is fairly obvious that under certain conditions we may 
exchange the role of the lengths L2 < Li in the derivation, 
such that the summation is carried out for the larger length 
Li instead. Thence the results (22) and (24) will remain valid 
upon swapping all occurrences of Li and Li if the aforemen- 
tioned conditions remain satisfied; we refer to this as the 
'swapped-lengths' version. As then, Li <C ^ is necessary, we 
are restricted to the small-lattice regime; the remaining con- 
dition is equivalent to Li/Lj <C 1, as follows from careful 
inspection of the B28 derivation of [11]. 

A consistency check between the result as it reads above 
and the swapped-lengths version has to take into account 
that it can only hold in the small-lattice regime, and that 
the orders of several terms change. Numerical comparison 
of both versions and comparisons with the evaluation of the 
exact result />rw show that, indeed, the swapped-lengths ver- 
sion reproduces the exact result more accurately than the 
form given above, as it circumvents the aforementioned con- 
vergence issue. This is only a minute problem however, and 
both versions agree and are in good agreement with the ex- 
act result (for the square lattice) over a certain range of dis- 
tortions (checked for ji = L2/L1 = 1/4 ... 1 for the small- 
lattice regime with 5 = 4-10^, a/W =10*, the agreement of 
the swapped-lengths version with the exact result extending 
farther down to ~ 1/10). 

5.4 Small distortion 

We now want to approach the question of the aspect ratio 
influence not from the extreme case, but in contrast starting 
to deform a quadratic lattice. Concretely, we ask when (and 
to what effect) the random walk feels the global change if 
we shghtly distort a quadratic lattice of 5 = sites, keep- 
ing W I a and S fixed. The latter condition is necessary since 
(judging from previous results) the effect of varying S is 
dominant otherwise; this hinders us from simply distorting 
the lattice lengths to "nearby" integer dimensions ( (L — 1 ) (L -|- 

Therefore we proceed as follows: We treat prw as a func- 
tion prw(M)5')§)» keep the latter two arguments fixed, and 
assume that the aspect ratio \i can be varied continuously. 
Any statements about the local behavior of that function 
around the quadratic shape value [i* remain sensible as long 
as higher-order terms do not restrict their validity to intervals 
too smaU. We do not believe this to be the case, basically 
because we cannot conceive of any mechanism that could 
render the behavior of pvti\^ ,^=const ii<ii* non-monotonic. 

The effect of a small distortion could depend on the regime 
we are in. We now focus on the small-lattice regime, since 
we expect aspect ratio effects to be most prominent there, 
and since for this case, we could obtain the general asymp- 
totic behavior in Section 5.3. Note that starting from a square 
lattice with the question at hand rules out the intermediate 



regime, and that of large lattices was argued not to be of 
much interest above. 

Now d V/ dAi ^ = [5( 1 - ^ ) ] - ' , and at the ex- 

tremal point, ^^qjP^ = _(p*(0;^V^^5^, since 
the first derivative vanishes (dp/d/i|^^^, = due to L\ <-> 
L2 symmetry of p). For our purposes we only need to find 
the sign of the second P* (0; ^ ) derivative. To this end, we 
use the approximation of Section 5.3, supposing that the 
dominant contributions to P* (0; ) are also dominating the 
sign of the second derivative. As a function of \l, (22) reads 



P*(0;,^) = const - 



M/3+5f/r 

r;r(l-2^o) 2(1-2^0) 



+ 



(3T)-l)7r 
36r ^ 



2(1-2^0)5 



+ 



^(L-^(l-^)l/2) 



(25) 



and therefore in the second derivative, only the ^3 contribu- 
tions and those from the square root term remain, viz. 



dV(0;^) ^ 



-5W/5)/dM^ 



dM^ 

^(1) 



2r(l-2^o) 



(26) 



The 5^' contribution is irrelevant due to the S ' factor, and 



2c(0) 



dAi^ 



167r (e-^'^" +6e-4'^'' +12e-^'^'' + . . . ) 



(27) 



where (as before), we only treat case (a) with r—\ = r\. Nu- 
merical evaluation of this term dX\i — \ yields approximately 
0.095. Thence, ^^^^^ > 0, and consequently, d^p/d/i^ < 
0: the encounter probability of the random walk in the small- 
lattice regime has a maximum for the quadratic lattice, and 
it decreases with any distortion from this shape.^ This is the 
quantitative foundation of our heuristic arguments in Sec- 
tion 5.6. 



5.5 The quasi-one-dimensional limit 

Finally we want to clarify the nature of the transition to ef- 
fective one-dimensionality: What happens if the aspect ratio 
goes to infinity, so that the macroscopic structure of the lat- 
tice becomes quasi-one-dimensional? 

First, we comment on the relationship between the re- 
gimes of the one-dimensional and the distorted two-dimensional 
lattices. In the 2d-large-lattice regime, £ is the smallest length 

^ We checked by numerical evaluation that P* (0; ^ ) is monotonic, 
or equivalently, that the first derivative of P*(0;^) indeed does not 
change its sign from ^ = 1 to /i = 1/4, corresponding to the distor- 
tion Li — > 2Li , L2 5L2 from a quadratic lattice with even lengths 
L = L\ = L2. Here we used the swapped-lengths version of (22) to 
avoid convergence issues and improve compliance with the exact re- 
sult as detailed in Section 5.3; the necessary conditions for this to be 
allowed are easily satisfied. 
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scale; this cannot be sensibly related to any one-dimensional 
situation. It is then clear that the large-lattice regime on a 
one-dimensional lattice of length £u corresponds to the 
intermediate regime of the two-dimensional situation, i.e., 
1 <C ^2 -C ^ -C Li, the common feature being that exactly 
one lattice dimension is much larger than the typical ran- 
dom walk. The remaining small-lattice regime of the two- 
dimensional case (1 <^ L2 Li £) has its analogon in the 
small-lattice regime 1 <C L <C ^ of the one-dimensional case, 
in both cases defined by a diffusion length that exceeds any 
lattice dimension. 

We now consider a scahng limit of the three lengths in 
the (two-dimensional) small-lattice regime so that we can 
apply the result of Section 5.3. Since £ is the largest length, 
we use ^ ^ °o as the basic scale in powers of which we ex- 
press the Li behavior. In the original form (22), 1-2/^1 ^0 
is not permitted, whence we employ the swapped-lengths 
version again. Now let L\ ^ I and L2 ~ i^l'^ (for the latter, 
any exponent in [1 /3 , 1 ) can be chosen, the lower bound pre- 
sumably an artefact of the details of the involved approxima- 
tions only). This implies L\/Li, ^ (r^l^ ^ 0, thus satisfying 
a central assumption for the approximations, while Li /L2 ~ 
^1/2 00, and consequently, the sf'^'^^'' terms vanish as well. 
Moreover, LiL2(l - ^ ^"'^^ ^ and L1/L2 •LiL2(l - 
1^ ) — > const. This scaling describes a situation where all lengths 
diverge, but the larger lattice dimension scales with the ran- 
dom walk length, while the smaller length increases more 
slowly to distort the lattice to vanishing aspect ratio L^/Li. 

Putting all of this together we obtain 



lnL2 



L1/L2 \ 



(28) 



where higher-order terms in the brackets were omitted. Here, 
the logarithmic correction terms (and coequal ^(1) terms in 
P* (0; t, )) can be seen to die out, while the second correction 
term actually approaches a constant. Neglecting all higher- 
order terms we thence have 



1- 



6(1-2^0)' 



(29) 



or for the natural choice of a square lattice (^0 = 1 /4), and 
with 1 — ^ w to leading order. 



(30) 



We compare this to the corresponding one-dimensional small- 
lattice regime result (62), which implies (cf . Appendix E) 



Prw,ld ' 



1 



(31) 



The proper scaling limit results in a crossover from loga- 
rithmic correction terms in the full two-dimensional small- 
lattice regime to the squared length ratio correction we found 



for the corresponding genuinely one-dimensional small-lattice 
regime. Only the larger length (scahng as the random walk 
length) remains, while the smaller one has disappeared from 
the result. Note that the appropriate rescaling of the random 
walk length, ll\p2. (accounting for the sphtting up 

of the number of steps between the two dimensions, as ex- 
plained in Section 5.2 and in Appendix E) reproduces the 
exact numerical pre-factor of the leading correction term, as 
it should. 

5.6 Discussion 

We now explain in a coherent fashion the main effects which 
govern the behavior of the encounter probability on a dis- 
torted lattice, as shown in Figures 4 and 5. We continue to 
keep S and ^ (or I) constant. 

The total decline from the peak value of prw to its mini- 
mum for the fully distorted S ■ 1 lattice strongly depends on 
the absolute lattice size: While for the absolutely large lat- 
tice, the probability dropped to less than 0.5% of that in the 
quadratic case peak, we still have 1 /6 of the peak value left 
when L2~ I with the absolutely small lattice. The quadratic 
lattice has either p ~ 1 or p ~ {a/W)/S, but upon maxi- 
mal distortion it invariably ends effectively one-dimensional 
(with halved hopping rate) in its 1 d-large-lattice regime, such 
that pid ~ \fafW /S, see Appendix E. Fixing the quadratic 
lattice regime via SW /a then impUes that p also stays fixed 



as a function of 5, and hence yields the ratio p/pid°^ S ^JWJa = 
y/SWjaVS VS in both regimes. With increasing length 
the two-dimensional walk misses many sites [24], but still 
sweeps an area of order i'^ in the large-lattice limit, while 
the one-dimensionahty of the extremely distorted case (also 
in its large-lattice regime) changes the £ power to the un- 
favorable. Obviously, a transition from a quadratic 2d-small- 
lattice into the Id-small-lattice regime would only yield minute 
corrections to near-perfect encounter probability. 

As for the shape of the prw dropoff, there are two effects 
which cooperate in decreasing the encounter probability due 
to distortion, related to the dynamics and the initial condi- 
tions of our problem. 

Regarding the dynamics, i.e., the exploration of the sur- 
face, consider the effect of 'wasted steps' (described in Sec- 
tion 5.2) in the smaller lattice dimension. Close to the quadratic 
shape there are no wasted steps, while close to maximal dis- 
tortion the effect is most pronounced, as steps in the smaller 
direction are nearly useless indeed - but even there the effect 
is that of halving the hopping rate, or modifying the random 
walk length £ hy a factor of ^(1). Expecting only power- 
law behavior of p in the lengths, this hardly shows on our 
logarithmic scale. Motion in the smaller dimension L2 does 
become 'wasteful' once £ ^ L2, when more steps or further 
reduction of L2 do not help the walker to get closer to the tar- 
get: Its residence probabiUty has already spread out in this 
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direction, and new territory can only be explored by step- 
ping out in the other dimension (Li). But in a regime where 
this substantially worsens the chance to reach the target, it 
is owed to extending the larger lattice length only. Hence 
the 'dynamic' effect lies in the changing ratios Li^i/^ upon 
distortion, governing how fast or to what extent the lattice 
is explored in each direction, and eventually switching the 
(refined) regimes. 

The second effect is 'static' in that it depends on the 
initial conditions and not on the dynamics of the system. 
We always assume deposition of the walker homogeneously 
distributed over the lattice. Even when, in the small-lattice 
regime, the largest parts of the lattice will be swept by most 
walks, it is still an obstacle for the individual walk if it starts 
further away from the target. Now as soon as the lattice be- 
comes distorted, the complete distribution of the walker's 
initial distance (and in particular its average) to the target is 
shifted to larger distances, as one can convince oneself of 
with a simple sketch. For the large-lattice regime, a region 
of an extension £ surrounding the target is not directly af- 
fected by distortion. But the spatial probability distribution 
of an immortal walker is then basically a Gaussian spread- 
ing with time, and relegating sites to a slightly further dis- 
tance greatly reduces the chance for a walker starting there 
to reach the target in a given number of steps. Hence rare 
long and successful walks are additionally suppressed by 
shifting the initial distance distribution. 

The relative importance of these effects is not easily quan- 
tified in general, as it depends on the absolute lattice size, 
the regimes, and the aspect ratio of the lattice, but we will 
explain their influence in the different regimes of Li 2 and 
£. Any numerical ^(1) factors in the comparisons will be 
omitted. 

First, let £ <C ii,2, the standard large-lattice regime in 
which the walk is always of two-dimensional nature. As ar- 
gued before, the different lengths of the lattice hardly have 
any effect here, and the encounter probability is governed by 
the ratio p ^ £^/S <^ 1 , independent of the aspect ratio . This 
region corresponds to the small peak plateaus at the right 
of Figure 5. In this regime, only the static effect might be 
important. 

For the small-lattice regime Li^2^£, the encounter prob- 
ability is close to unity. A small fraction of walk reaUzations 
does not lead to recombination, but in this regime, avoid- 
ing the target does not become much easier when the lattice 
is distorted. Now the aspect ratio determines whether the 
random walk behaves essentially one- or two-dimensional. 
Based on (28) of Section 5.5 we expect the crossover scale 
between the two types at an aspect ratio of roughly n = 
l/lnL2 below which (small-lattice) one-dimensional behav- 
ior prevails, namely /? « 1 — (Li /£)^/3. The main difference 
to the two-dimensional behavior (apart from the absence of 
logarithmic corrections) is that only the larger length still en- 



ters. But the respective term has to be small anyway for the 
expression to be applicable, and we cannot expect to see it in 
the plots. We thus only have a minute effect of the distortion 
in this regime as well, a fact easily seen at the peak plateaus 
in Figure 4. Actually, the plateaus are more pronounced in 
this case than for the corresponding large-lattice plots, be- 
cause not even the static effect has any significant influence 
anymore - basically the whole lattice is swept anyway. 

Finally consider the intermediate case L2 ^ f < Li . This 
regime is more complicated since we do not know whether 
dying or meeting the target is dominant in setting the resi- 
dence time of the walker and the magnitude of p. The as- 
pect ratio determines whether we deal with a genuinely two- 
dimensional lattice, or rather with a lattice so elongated that 
it is of one-dimensional nature. 

If jU = L2/L1 ^ 1 is extremely small, we can effectively 
consider the system as one-dimensional and in its Id-large- 
lattice regime. Homogenization in the L2-direction is much 
faster than the spreading in the Li -direction then, and thence 
the probability to reach the target is essentially that of reach- 
ing the projection of the target position onto the L\ -dimension. 
This roughly coincides with the probability for the walker to 
start within a reach £ of the projected position, and this is 
p ^ £/Li. For the effectively one-dimensional intermediate 
lattice this reasoning is justified, as is shown by Appendix E. 

In contrast, if the aspect ratio is not smaU enough for 
this viewpoint, the nature of the two-dimensional random 
walk shows: it spreads without fully exploring the swept 
area (rather with increasing 'sponginess' of the set of vis- 
ited sites, again a specialty of spatial dimension two). Con- 
sequently, we rather expect anS/£^ dependence with typical 
logarithmic corrections as for the quadratic case. L2/L1 is 
only moderately small now, so we may use the result of Sec- 
tion 5.3 provided that {Li/ £){Li/ £) <c 1. From this we get 
p « 1 — (5/£^)ci (const. -|-lnLi), a functional dependence 
similar to the small-grain expansion for the quadratic lattice 
case. The crossover aspect ratio between both types of be- 
havior could not be determined, because in this regime, the 
expansion breaks down on the way to the one-dimensional 
asymptotics. For the quadratic periodic lattice, the logarith- 
mic correction emerges from integration of the slowly decay- 
ing return probability of long walks; with a distorted lattice, 
the dominant bounding contribution stems from the larger 
length only, as is most easily seen by rescaling the approxi- 
mate integral expression. Appendix C. 

This intermediate regime governs the behavior of the 
Figures 4 and 5 once we leave the plateau around the quadratic 
shape. In all cases, we first enter a linear decline on the 
log/i scale, as predicted for the two-dimensional intermedi- 
ate regime (^ — InLj ^ log jtt with S constant). This behav- 
ior roughly starts once the larger (smaller) length of Li_2 be- 
comes of the order of I, depending on whether we started in 
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the small-(large-)lattice regime, which is precisely the con- 
dition to enter the intermediate regime. 

For the leftmost (most distorted) part of the plots, the 
absolute size of the lattice becomes crucial. On absolutely 
large lattices (thick lines in Figures 4 and 5) the linear de- 
cline ends in an exponential shape close to the fully distorted 
lattice. This is explained by the effective one-dimensionality 
as argued above; we end up with the Id-large-lattice behav- 
ior p r-j The crossover aspect ratio is read off to be 
roughly between 1 /400 and 1 /lOOO. 

In stark contrast, the lattices of small absolute size (thin 
lines in Figures 4 and 5) do not show any clear deviation 
from the 2d-intermediate-regime decUne down to full distor- 
tion. First, the total decline is not as extreme, as explained 
before. Second, we still are in the Id-limit for full distor- 
tion. But leaving ii = l/S, immediately two-dimensional ef- 
fects and the corresponding logarithmic terms (linear on our 
scale) dominate the p behavior, because the aspect ratio is 
far less extreme than for the absolutely large lattices - con- 
sistently, the crossover aspect ratio determined there is not 
reached here. 

6 Conclusions 

In this work, we have thoroughly examined the encounter 
probabihty of two mortal random walkers on a periodic lat- 
tice as a model of a confined geometry. We compared the 
results of continuum diffusion models, exact random walk 
treatments and their asymptotic behavior and heuristic re- 
sults obtained from standard random walk analysis. We high- 
Ughted their similarities as well as explained the crucial dif- 
ferences, most importantly the features responsible for the 
slow logarithmic convergence to the continuum Umit, which 
is a direct consequence of the criticality of dimension two 
for diffusion and random walks. The discrete-time random 
walk results have been shown to be in excellent agreement 
with kinetic Monte Carlo simulations of the continuous-time 
version that naturally lends itself to physical applications. 
On a side note we could corroborate earlier claims that the 
lattice type used is of minor importance to all results dis- 
cussed herein. 

The second half of this article has been devoted to the 
analysis of the influence of the geometry of the lattice, ex- 
amined at the example of an aspect ratio differing from unity. 
We considered an extremely distorted lattice and explained 
our findings for this situation. Then we generalized an early 
result by Montroll to the distorted situation, and thus could 
examine the effect of distortion starting from a quadratic lat- 
tice. Moreover, we determined a scaling limit in which the 
dying out of logarithmic terms (characteristic for the two- 
dimensional situation) in favor of algebraic corrections (typ- 
ical of the one-dimensional case) can be nicely seen and in- 
deed, the one-dimensional asymptotic behavior is fuUy re- 



covered. Finally, we gave a general explanation of the ef- 
fects which govern the changing encounter probability and 
the regimes to which they pertain. 

We believe that this work fills a gap in the vast litera- 
ture on random walks and first-passage problems. Most im- 
portantly, we have achieved the goal stated in the abstract: 
We have shown that the (completely analytic) treatment of 
the simple discrete-time random walk model provides a full 
quantitative understanding of the encounter probability for 
the (homogeneous) continuous-time reaction-diffusion sys- 
tem as simulated by the kinetic Monte Carlo method, and 
this holds for all the aforementioned aspects including the 
shape of the lattice and its structure. These results are impor- 
tant for many applications ranging from astrochemistry to 
biophysical problems, and they equally matter for both the 
fundamental theory as well as for the numerical simulation 
of reaction-diffusion systems. 

Further work will be centered on the role of quenched 
and annealed disorder in the rates of hopping and dying of 
the walkers, mainly from the point of view of the recombi- 
nation efficiency of systems such as those described herein. 
An important next step will be the assessment of the validity 
of the master equation framework in the homogeneous sit- 
uation, where hitherto, spatial correlations between random 
walkers on the lattice have been tacitly neglected. 
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A Derivation of prv, 

We closely follow and mimic the notation of [17] without citing indi- 
vidual known results. All of the techniques used here have been devised 
quite some time ago, see e.g. [16] with the exact same formula as (42), 
and we basically put together all the necessary pieces in an appropriate 
form. 

We deal with a finite periodic homogeneous lattice with S = YY'i= \ 
sites in total, or a rf-dimensional torus with extensions Lj in the 7th di- 
rection. Onto this lattice, we put two random walkers (in discrete time), 
starting at random sites, independently and homogeneously distributed. 
They do not interact except when they meet. The walkers are assumed 
mortal (corresponding to desorption) with constant and equal survival 
probability S, per step. Our question is: "What is the probability that 
these two walkers meet before one of them dies?" 

This problem can be mapped to that of a single walker, starting 
from a random site sq, with the same survival probability per step, the 
question being with what probability it eventually reaches a certain 
fixed site s* on the lattice without dying prematurely (all sites are la- 
beled by a variable 5 whose structure is irrelevant right now). 



Dying and moving of the walker are independent. Thence 

PrjMortal walker reaches site s* for the 
first time on the feth step} 
=Pr{Mortal walker has completed 

at least k steps} x (32) 
X Pr{ Immortal walker reaches s* for the 
first time on the fcth step} 

= <^*-Fi(.s-*|,s-o), 

where Fi^{s\sq) is the probability (of an immortal random walker) of 
arriving at site s for the first time on the k\\\ step, given that the walk 
started at site .so- It should be noted that we adopt the convention that 
^o('5|'So) = 0> and therefore do not count a walker already starting at .so. 
We do not care for the time when this first passage of s* occurs, and 
are therefore interested in the quantity 



F{s*\so;^):=Y^^%(s*\sa), 

k=0 



(33) 



which happens to be the generating function of Fi^{s*\so). One can con- 
vince oneself that, by the definition of as the ^raf-passage proba- 
bility, every encounter is counted only once, and so F(.s*|.so;i^) is the 
probability of a mortal random walker with survival rate t,, starting 
from .So, to reach s* before dying. 

For the remainder of the derivation we only have to be concerned 
with immortal (usual) random walkers. Then for any random walk, 
there is the relation (e.g. [17] Eq. (3.27)) 



P{sW,E,)-8s, 
P{s\s;^) 



(34) 



where P{s\sq;^) is the generating function of /^fc(i|io)> the probability 
that the random walker, starting at site ^q, is at site i on the fcth step, 
with the convention that PoC-^ko) = 5s,s„. 

However, regarding the convention used for the first-passage prob- 
ability, we want to count a walker that starts at the target site as "reach- 
ing it for the first time" on the 0th step, and not count later returns. 
Denoting the corresponding probability by Fii{s\sQ), it is related to the 
original one by 



hiA^o) = Ss,soSk,Q + (1 - Ss,so)Fkis\so) ■ 
For the generating functions this implies 
F{s\so;i) = Ss,so + (1 - 5s,jF(s\so;^) , 
and inserting (34) yields 



P(,s|.s;^) 



(35) 



(36) 



(37) 



Now we can already write down the answer to our problem in very 
general terms. What we actually want is an average S^^ Y.soen PiA^o) 
of the revised first-passage probability (37) over all starting sites Sq in 
our finite periodic lattice Q, giving us the encounter probability of the 
original two random walkers: 



SP(s\s;^) 



(38) 



The apparent dependence on the target site .s will vanish in passing to a 
homogeneous setting. 

For further evaluation, we need to obtain P(.s|.so;i^). Let us for 
the time being denote by this quantity the generating function of the 
residence probabiUty of a random walk on an infinite lattice. Further- 
more, otir walk is homogeneous or translationally invariant, and we 
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may therefore use a single vector 1 pointing from the starting site .vo to 
the final site s as our variable. Such 'lattice vectors' have d components, 
all of which take arbitrary integer values that can be thought of as com- 
ponents of a 'spatial vector' with respect to the d fundamental vectors 
of the lattice. Then it is well-known that (e.g. [17] section 3.3.1) 



1 



I 



exp(— ;lk) 
l-^^(k) 



(39) 



with the integration domain the first Brillouin zone. A(k) is the 
structure function of the walk, defined as 



A(k)=X;exp((Tk)?{l), 



(40) 



the sum running over all lattice vectors, and q{\) being the probability 
of a step translating by 1. Since = 1, |'^(k)| < 1. 

The transition to a truly finite periodic lattice is now easy: Attach- 
ing a star to the respective residence probabilities, we have 



^;(')=L^n(l + Lm) 



(41) 



where L = diag(Li , . . .,Ld). The sum runs over all translation vectors 
of the (infinite) lattice, and this implies a completely analogous relation 
for the generating functions. Obviously, P* (1 + Lm) = P* (1) for an ar- 
bitrary m in the infinite lattice. From now on the vector 1 is understood 
to lie in the subset Q of the infinite lattice that stands for the finite peri- 
odic lattice (we use the same symbol for both the finite sets of sites and 
of associated translation vectors). It can be shown that consequently 
we have^ 



1 ^ exp[-2;i:(l(L"'iii) 



(42) 



The sum over the finite lattice Q is explicitly a multiple sum over the 
components of the vector m, each ranging between my = . . .Ly — 1. 

Rewriting (38) for the homogeneous walk and substituting the ap- 
propriate starred probabilities, we obtain 



SP*(0;^) 



(43) 



The numerator of this expression is nothing but the ^ transform of 
Elefl Pn (')' but this is unity due to conservation of the immortal walker, 
so that 



le£2 



1 



1-4 



(44) 



Again using (42) for the denominator of (43) jaelds a general expres- 
sion for the encounter probability on a homogeneous regular finite pe- 
riodic lattice, namely (1) of the main text. 

It is often desirable to consider a different convention that does 
not allow both walkers to start at the same site. To this end, one simply 
has to exclude the lattice distance (when talking about translationally 
invariant walks) from the average (38), the result of which is given in 
the main text as well. Note that this is not equivalent to simply employ- 
ing the original first-passage probability F instead of F, which would 
allow this situation, but not appreciate it as an encounter, and would 
simply lead to ^ instead of p„. 

^ The underlying identity reads £„exp[— ;(Lm)k]/(2;r)'' = 
— 2;rL~^m)/(n;Ly) and works component-wise. In a way, 
this is a more general expression than that for the infinite lattice, which 
can be recovered by sending Lj — > <x>. 



B Evaluation of one sum 

We now restrict ourselves to a still fairly large class of walks, namely 
those with structure functions 

A(k) = 2q(i[coski -|-cos^2] 

+ 2{cj\ +c/2)cos/:i cos^2 +2(172 — <?!) sin/:i sinfo ■ (45) 

Such a structure function belongs to walks with transition probabilities 
qo to take a step into either lattice unit direction, and q\ and q2 the 
probabilities to step into direction (1,1) or (—1,-1), and (1,-1) or 
(—1,1), respectively, subject to the normalization 4qo+2(qi +q2) = 1. 
Clearly, this includes the aforementioned cases of the main text: The 
isotropic square lattice 'type (a)' corresponds to qo = 1/^ and qi = 
(?2 = 0, and the isotropic triangular lattice walk 'type (b)' is represented 
by <?o = g2 = 1/6, qi =0. 

For this class of walks, one summation can be explicitly evaluated 
[11]. The result can be generalized to the Li ^ L2 case by some easy 
accounting work and then reads 



P*(0;4) = 7^- L [l-2^o<^cos(2;rmi/Li)]-' 

^1 mi=0 

1 2L2 

X [I_p2j-V2^ , 

1 2x^1 COsZ.20m[ "l~-^mi 

Here, < x„, = [1 - (1 - p,?,,)' ^^J/p,,,, < 1, and with 
2^ [qo + {qi + qi) cos{27lmi/L[ )] 



(46) 



W2 = 



1 — 2<^^ocos(2;rmi /Li ) 
2^{qi -^2)sin(2;rmi/Li) 
1 — 2i^^ocos(2;rmi/Li) ' 



(47) 



< < 1 and (j>„^ are given by wi -|- iw2 = pmi e"^"^ ■ In general, this 
yields 

Pmi =wi+wl 

(2',!- [</0 i '/T i I 2(/|(/2a)s-^^ ^ 2c/„lc/, I (/2 ) cos -^-i-J 



[1 -2^5ocos(2;rmi/Li)]2 



and 



(51 -92) sin 25a 



90 + (?i+92)cos2^ 



(48) 



(49) 



(whenever well-defined). 

For case (a) one obtains p„^ = [2/^ — cos(2;rmi/Li)]^' and = 



0, such that 

For type (b) one has 

v^y I I cos(2;i:;;i|/L| ) 2|Cos(7r;«i /L| )| 



3/^ -cos(2wmi/Li) -cos(2;rmi/Li) 



(50) 



(51) 



and 



sin(2;rmi/Li) 

= - [l+cosi2nm,/L,)] = (52) 

(whenever well-defined), which yields = —7l{mi/L\ — \2in\/L\\). 
The peculiar form of the angle is necessary to assure \^„^ | < 7r/2 cor- 
responding to Wl > 0, and while the chosen expression may result in 
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a wrong sign of sin^,,,, , only the unaffected cos{L2^mi) is used in the 
remainder. Hence we obtain 



100 -Sp/p 



P*(0;^): 



3/^ 



Li4^o2|cos(;r/ni/Li)| 



2L2 



(53) 



1 - 2x^1 COs(L2(|)mi ) ' 

and throughout, x„j (p„j ) as given above. 

The alert reader may object that, in case (b), there is actually one 
term for which both w\ = = W2 and thus p,,,, vanishes (see our valid 
explicit result for the latter), namely for even Li and mi =L\/2. Hence, 
x„| , <|),„| and the last expression for P* (0; ^ ) are ill-defined then. From 
the definition of jc,„| one can see that for — > 0+, r; p,„^ /2 ^ 0+ 
as well, and the second line of (46) (corresponding to the evaluated 
'inner sum' in the derivation of [11]) converges to unity, which is the 
correct value of the original quantity. Our last expression complies with 
this via canceling singularities - we chose the simplest form of the 
result, which has to be sUghtly altered for numerical evaluation. 



C Large lattice approximation 

For all approximations in this and the following Section it is justified to 
treat 1 — ^ and (1 — ^)/^ = W/a synonymously due to ^ < 1, which 
will no longer be mentioned when it only introduces higher-order er- 
rors compared to the desired accuracy. Moreover, we will treat loga- 
rithms as being of the order of unity. This might make some expres- 
sions more cumbersome, but it is numerically adequate. 

By letting Li^2 — > ""inP* (0;<^) of (1) one obtains a double integral, 
which after two linear substitutions reads 



(1-^)5 



1 



(2;r)2 



/ 

[0,2jt]2 



dMl dU2 

l-^A(u) 



(54) 



X periodicity once again allows us to shift the patch of integration to 
the Brillouin zone B defined in Appendix A. As is shown in the 
Appendices of [17], in terms of the complete elliptic integral of the 
first kind. 



1 

K{k) = [ ^ , with \k\ < 1, 

one can derive that P(0;(^) = ^K{^) for the square, and 



X K 



2{c+- 



;c- + i)(c+-i) 



(55) 



(56) 



for the triangular lattice, where c± = 3/(^ + 1 ± ^/S + G/E,. Now we 
need the expansion of the elliptic integral for k < i (which is also the 
case for the triangular lattice if ^ itself is close to unity), viz. 



n=0 



2nk 



-(l/2)ln(l - k^) + Yin+\) - Vf(n+ 1/2)] 

Fa)/r(a) is the POCHF 
a{a -I- 1) ... (a -I- n — 1) for positive integer i 
in the exact results for P(0; ^ ), one obtains 



, . (57) 

here (a)„ = r{n + a)/r(a) is the Pochhammer symbol, equal to 
(a-l-n — 1) for positive integer n. Using these expressions 

iilfc fnr P( n- ^ nnp nhtninc 



P(0;§) = [l + ^(l-^)]x 



iM8/(l 



^ln[12/(l 



This yields the expressions for 



square lattice, 
] triangular lattice. 

(58) 

presented in the main text. 




logs 



Fig. 6 Relative difference of pmc w.r.t. for SW /(4a) = 10 type 
(a) (dashed / squares), type (b) (dotted /triangles). 
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Fig. 7 SW/i4a) = 10"^ 
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Fig. 8 SW/(4a) = 10"^. 
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D Simulation results for 

Figures 6-12 show the relative error (per cent) 1 00 • [pmc — Prw ) / Prw as 
a function of the lattice size S, for quadratic lattices. Lattice sizes and 
rate ratios a/W are chosen as described in Section 3.3, so that one plot 
roughly corresponds to a fixed regime and about constant p, to keep the 
standard deviation of the same order of magnitude. Outside these pa- 
rameter ranges, nothing interesting happens; in the large lattice regime, 
the leftmost data points are omitted as they no longer satisfy a/W » 1 . 
We have plotted a corridor of half width G around perfect coincidence 
to show that the discrepancy between analysis and simulations is statis- 
tically insignificant. 

Figures 13-16 show corresponding simulation results for rectangu- 
lar lattices of varying aspect ratio. Here S is constant for one plot. We 
restrict ourselves to S = 400 and S = 4 x 10* as examples of absolutely 
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Fig. 10 SW/(4a) = lO". 
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Fig. 11 5W/(4a) = 10'. 
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Fig. 14 S = 4 X 10^, W/a= W^. 
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Fig. 15 S = 4 X 10^ W/a = 10"^ 
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Fig. 16 S = 4 X 10^ W/a = 10"^ 

E Asymptotics of the truly one-dimensional result 

For the one-dimensional symmetric random walk with structure func- 
tion X{k) = cask and a lattice of 5 = L sites, we have 



no;^)=-i:^ 



1 



Q 1 — cos{27lm/L) 



(59) 



The identity of Appendix A in [1 1] (with wi = , 1^2 = 0, so that p„j 
^ and (j)„i = 0, and ;c„j = [1 - (1 - ^^Y^^]/^) then yields 



^L+[l_(l_^2)l/2]i 



Fig. 13 Relative difference of pn,(. w.r.t. pr„ for 5 = 4 X 10 , H' /a = V"'S; _ ^2 ^\ — ^i- — [1 — (1 — ,^2)l/2]t 



10 , type (a) (dashed/ squares), type (b) (dotted / triangles). 



small and large lattices, respectively. Further, we keep the "regime" (in 
the original sense introduced for quadratic lattices) constant as well: 
for S = 400, we choose a/W = 100 and a/W = for S = 4 x lO*" 
we test with a/W = 10* and a/W = 10^ (the latter figures belonging to 
the small-lattice regime), but note that the refined regimes for Li=^Li 
actually change when distorting the lattice. We have chosen our val- 
ues so that the larger length L\ exceeds the random walk length (. in 
the large-lattice regime throughout, and that it is much smaller than 
this length in the small-lattice regime for the quadratic case, but finally 
becoming much larger than it when distorting the lattice. Again, we 
show the relative error of simulations with respect to the random walk 
result in per cent, and plot a standard deviation (as obtained from prw) 
corridor for comparison. 



(60) 

For comparison with the asymptotic two-dimensional behavior, we 
evaluate the asymptotics of the Id-expression for a 'small' lattice, i.e. 

for L < ^ = v'(^/(l -t,). Let a = y'l - t,^ < 1, and note that this is 
not exactly the a of [ 1 1] . However, L/£ r; L\/ 1 - t, = L\/ 1 - \/l - f 
La/ \/2, and since this is to be ^ 1, we can still use La ^ 1 in the fol- 
lowing. With 1-^ = 1- VT^^ and pj = L{1 - ^)P*{0;^) we 
then obtain after minor manipulations 

, L[l-(l-«^)'/^]l + (|^)"^^ 



a 1 _ /i-oi\V2 

^ \l+a) 



1+0 

Straight-forward expansion in La ^ 1 yields 

Prx, Id = ^ "I 



-1 -i+(^ + ^(La)3, 



(62) 
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where we omitted all terms of relative order a. This is sufficient for 
the sought scaling limit L ^ oo, a ^ Q while La ~ const., and it is 
additionally justified by a Mathematica check of our calculations, 
which shows that the only interesting a orders left out in the result are 
terms of ^(a^) and thus much smaller than {Ldf = L^(l - « 



Obviously, the most important difference compared to the 2d case 
is the absence of logarithmic terms, which are a characteristic sign of 
the marginal dimension two. 

For completeness, let us also give the large-lattice asymptotics, i.e., 
for the case La 2> 1. Starting from the still exact (61), we use (1 — 
«)/{! + a) = 1 - 2a + ff(a^). For La > 1, this raised to the L/2 
power is much smaller than unity, the correct expansion thus reading 
exp(— La)(l + ff{La^)), and here we assume that terms of the latter 
relative order can be omitted as being much smaller than unity - this 
is an example that may be refined as necessary. One thus obtains the 
large-lattice restilt 



Note that in this regime, the ratio of lengths no longer appears squared. 
In one dimension and for large lattices, the essential question for the 
encounter probability is whether the walker is deposited in the range t 
from the target, and the probability for this to happen is the length ratio 
raised to the lattice dimension. The difference to the two-dimensional 
case is that the Id random walk explores a dense region instead of a 
sponge-like structure. 

Lastly, corresponding asymptotics for the extremely distorted ver- 
sion of the originally two-dimensional random walk are obtained by a 
mere rescaling t — > l/V^, as can be seen from Section 5.2. 



l{Llt)\ 




(63) 



